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We report on density functional total energy calculations of the step formation and interaction 
energies for vicinal TaC (001) surfaces. Our calculations show that double and triple- height steps 
are favored over single-height steps for a given vicinal orientation, which is in agreement with recent 
experimental observations. We provide a description of steps in terms of atomic displacements and 
charge localization and predict an experimentally observable rumpled structure of the step-edges, 
where the Ta atoms undergo larger displacements compared to the C atoms. 
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The stepped surfaces of transition metal carbides (TiC, 
TaC, etc.) are technologically important due to their ap- 
plications as surface catalysts with high corrosion resis- 
tance. From a fundamental point of view, these surfaces 
have received attention as they provide a rich ground 
for studying faceting transitions that involve single and 
multiple-height steps. In a series of experiments, Zuo and 
coworkers have investigated the structure of sev- 

eral surfaces vicinal to TaC (001) and found that multiple- 
height steps are favored over single-height steps in all the 
cases. Motivated by these experiments, we perform ab- 
initio density functional calculations to obtain the first 
quantitative results for the step formation and interac- 
tion energies for steps of different heights on these sur- 
faces. 

The preference for single or multiple-height steps for a 
certain vicinal orientation is determined by the energy re- 
quired to create an isolated step of a certain height and 
the repulsive interactions between these steps. In gen- 
eral, it costs more energy to create multiple-height steps, 
as the number of bonds that are broken in the process 
of creating a step increases with its height. On the other 
hand, the energetic contribution from the repulsive step- 
step interactions can be the lower in the case of a vicinal 
surface with multiple-height steps due to a lowering of 
the number of step pairs per unit length. As we will 
see, both of these effects are borne out by our calcula- 
tions. Furthermore, we will show that the increase in the 
formation energies of double and triple-height steps over 
the single-height steps is not large enough to offset the 
effect of step-step interactions. Consequently, the for- 
mer types of steps are favored for the vicinal orientations 
that we consider. In what follows, we also will report 
on the structures of different types of steps in terms of 
charge density and atomic displacements, and develop a 
phenomenological model for the electrostatic step inter- 
actions. 

The density functional calculations were carried 
out with the VASP package ||, using the projector 
augmented- wave pseudopotentials |J and the Perdew- 
Wang functional form for exchange correlation en- 
ergy. Bulk calculations with an 8-atom cell and a 35 
k-point Brillouin zone sampling gave an optimum Ta-C 
separation of a = 2.2415A. For surface calculations, the 
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TABLE I: Rumpling relaxation of TaC (001) compared with 
previous results ^|, Jiofc . The displacements of the first 
(dc,Ta) an d second {d 2 c ^ Ta ) layer surface atoms, and the rum- 
pling amplitudes n, r 2 (defined in p|) are given in A. 



Brillouin zone was sampled using a T-centered 8x8x1 
grid. This sampling yielded 15 k-points for TaC (001) 
and 21 to 25 k-points for the vicinal surfaces. The ions 
were relaxed using a conjugate-gradient algorithm until 
the total energy converged to O.OleV. The energy cut-off 
for the plane waves was set to 400eV (29.40 Ry) in all 
the calculations. 

Before we proceed with the discussion of vicinal sur- 
faces, we compare the predictions of the surface structure 
of TaC (001) obtained using the present model with ex- 
perimental observations and previous theoretical studies. 
It is well-known |^|-||lo[ that this surface exhibits a rum- 
pling reconstruction, where the C atoms are displaced out 
of the surface while the Ta atoms are pulled inwards. It 
can be seen from Table Q that our results for the relative 
displacements of the atoms in the first and second surface 
layers are in close agreement with available experimental 
H and theoretical |l(J data. 

For vicinal surfaces, we used super cells as illustrated in 
Fig. [l|, where L x and L z denote the terrace-width of the 
stepped surface and height of the supercell, respectively. 
The dimension of the supercell in the ^-direction (L y ) 
is determined by the period in the direction parallel to 
the step, which in the present case is the lattice constant 
of TaC, 2a. To create the steps, we employ shifted (i.e. 
non-orthogonal) boundary conditions jn|, in which the 
amount of shift in the z- direction determines the height 
of the steps. The structures of the single (Sh), double 
(Dh) and triple (Th) -height steps that we consider in 
the present work are shown in Fig. [l| All calculations on 
vicinal surfaces were carried out using a slab of TaC with 
8 atomic rows (17.932 A) and a vacuum thickness of 10A. 
In order to ensure that these results are not influenced by 
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FIG. 1: Geometry of a typical stepped surface, with the su- 
percell indicated by dashed lines. The structures of the single 
(Sh), double (Dh) and triple (Th) -height steps considered in 
the present work are also given in the figure. 
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268.89 998.42 119.96 


88.85 
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466.47 550.34 104.05 


6.11 
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637.08 1217.92 94.74 
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TABLE II: Step formation energies A;£ and repulsion strengths 
A^ obtained using linear fits of data in Fig. to Eq. (FD). The 
table also shows the scaled step formation energies and 
interaction strengths (3% introduced in Eq. (ph. 



finite-size effects, we repeated the calculations for certain 
selected step structures and terrace widths by increasing 
the cell thickness to 10 atomic rows and the dimension 
of the vacuum to 15 A; the results were found to change 
by less than 5% in all the cases that were considered. 

The energetics of stepped surfaces can be understood 
by considering the ledge energy, defined as A n = (E — 
N p eb — joA)/2Ly [ p] , where E is the total energy of the 
N p Ta-C pairs in the slab, — — 22.1536eV is the bulk 
energy per pair, 70 = 95.53meV/A 2 is the surface en- 
ergy of the TaC(OOl) surface and the subscript n denotes 
the step- height. The ledge energy is the energy per unit 
length of the vicinal surface in excess of the surface en- 
ergy of the terraces that separate the steps, and includes 
both the step formation and interaction energies. Using 
the fact that both the elastic and electrostatic effects give 
rise to dipolar interactions jl2|, Q , the ledge energy can 
be expressed as 
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where A n and A^ denote the formation energy and the 
dipolar interaction strength of a step with height na, re- 
spectively. The computed ledge energies for the three 
types of steps are plotted in Fig. ^ along with linear 
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FIG. 2: Ledge energies A for the single (Sh), double (Dh) and 
triple (Th) height steps plotted as a function of The 
points for which the step separations exceed a —2. 24 A are fit 
to Eqn. (0), A n = A£ + A^/L 2 Q. Parameters A£, A^ are 
given in Table [n|. 
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TABLE III: Surface energies (in meV/A 2 ) of the (102) and 
(103) surfaces considered by Zuo and coworkers in The 
table gives surface energies calculated using Eq. ^ and the 
scaled fitting parameters in Table |, as well as surface en- 
ergies for (102)-Sh and Dh and (103)-Sh calculated directly 
from the ledge energies marked in Fig. ^. The close agreement 
between the surface energies obtained using the two different 
approaches indicates that Eq. (Q) gives an accurate descrip- 
tion of the ledge energies. 



fits to the functional form assumed in Eq. ([l]) [[Til . The 
fitting parameters A n and A^ are given in Table ||. 

It can be seen from Table |l| that the formation energies 
of the steps increase with their heights, in agreement with 
bond counting arguments presented earlier. The dipolar 
interactions, on the other hand, do not show a monotonic 
trend; the interaction strength of the Dh step is smaller 
than the corresponding values for the Sh and Th steps. 
The relative importance of these two contributions can be 
understood by considering the surface energy of a vicinal 
surface made of steps with height na, given by 



7n(0)=7ocos0 + /#|sin0|+/# 



sin 3 1 
cos 2 
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where 6 is the vicinal angle and the parameters (3^ = 
A^/(na) and /3% = A^/(na) 3 denote scaled step forma- 
tion and interaction contributions, respectively. Using 
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the values of these parameters tabulated in Table |T[ we 
find that both the scaled formation and interaction parts 
show a monotonic decrease with increasing step height, 
indicating that multiple-step heights are favored over the 
single height steps for any vicinal orientation. In the fol- 
lowing paragraph we consider experimental work of Zuo 
and coworkers and show that the predictions of our 

calculations are in agreement with their observations. 

Using a statistical analysis (H| of the step separation 
distributions on TaC(103) and (102) surfaces, Zuo and 
coworkers have determined that the former surface con- 
sists of a large number of triple-height steps and a smaller 
number of double-height steps while the latter surface 
consists of almost equal amounts of areas that contain 
either triple or quadruple- height steps. In order to un- 
derstand these results, in Table III, we have used Eq. (||) 
and the data in Table || to calculate the surface ener- 
gies of these vicinal surfaces made up of Sh, Dh and Th 
steps. For both the surfaces considered in the experi- 
mental work, we find that the single-step vicinal is much 
larger (by ~ 10-20 meV/A 2 ) in energy than the vici- 
nals with multiple- height steps. Furthermore, the triple- 
height steps are favored over the double-height steps (by 
~ 4 meV/A 2 ) in both the cases, explaining their stabil- 
ity and their presence on (103) and (102) surfaces. The 
presence of a few Dh steps and the complete absence of 
Sh steps on (103) surfaces is consistent with the fact that 
the former is closer in energy to the Th steps than the 
latter by about 8 meV/A 2 (refer to Table At the 

present time, we are unable to consider the energetics of 
quadruple- height steps due to the size of this model [[l6| . 

We now turn our attention to the description of the 
steps in terms of the atomic displacements at the step- 
edges, which are shown in Fig. || for the (105)-Sh and 
the (102)-Dh surfaces. The figure shows that the atoms 
at the step-edges undergo both horizontal (x) and verti- 
cal (z) displacements. On moving a few atomic spacings 
away from the step-edge, the atoms on the terraces are 
predominantly displaced in the vertical direction, as in 
the case of the (001) surface. Overall, the displacement 
field of the step resembles the field produced by a force- 
dipole jl2| located at the step-edge, though there are 
new features that emerge. In particular, the step-edge is 
rumpled, such that the horizontal displacements of the 
Ta atoms at the step-edges are significantly larger than 
the corresponding displacements of the C atoms. Since 
both the components of the rumpling vector s given in 
Fig. H are comparable to the rumpling amplitudes for 
Ta(001) given in Table |, we believe that this effect could 
be observed in experiments. The step rumpling also leads 
to formation of electric dipoles, which will be considered 
next. 

The magnitudes of the electric dipoles at the step-edges 
due to rumpling can be estimated as Q|s|, where Q is the 
amount of charge transferred from C to Ta due to the 
partially ionic character of the TaC-bond. Since Q ~ |e| 
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FIG. 3: Atomic displacements of the TaC(105) and TaC(102) 
surfaces with Sh steps (top) and Dh steps (bottom), respec- 
tively. All vectors lie in the xz plane and have been magnified 
for clarity. Each atomic row is characterized by a rumpling 
vector, whose components are given by the difference in dis- 
placements of the Ta(red) and C(blue) atoms. The rumpling 
vector s of the atomic row at the step-edge is marked in the 
figure. The components of the rumpling vectors (in A) for 
a few atomic rows near the step-edge are also given in the 
figure. 
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FIG. 4: Electron localization function (arbitrary units) in the 
(010) cross-section for TaC(105)-Sh and TaC(102)-Dh sur- 
faces. The dipoles of the C atoms at the step-edges are de- 
noted by p. In the case of the Dh step, p' denotes the dipole 
of the C atom located below the the Ta atom at the step-edge 
in the plane immediately behind the one shown in the figure. 



and |s| ~ 0.1a (refer to Fig. g), where e is the electronic 
charge, the linear dipole density due to step rumpling is 
P ~ 0.1|e|. In addition to the dipoles that arise as a 
result of step-rumpling, we also find that there are ad- 
ditional dipoles that are created due to charge distor- 
tions at the step-edge. We observe that the valence shell 



4 



around the C atoms expands out of the surface leading 
to formation of electric dipoles, denoted by p and p' in 
Fig. ||. Such electronic distortions also exist for the Ta 
surface atoms, but to a much smaller extent. The mag- 
nitude of these dipoles can be obtained from the elec- 
tron localization functions, which are shown in Fig. ^ for 
the (105)-Sh and the (102)-Dh surfaces. Our calculations 
show that the typical magnitude the linear dipole density 
due to charge spill-out at the step-edge isP~0.3 — 0.7|e|. 

The interaction energy between steps characterized by 
the effective electrostatic dipole density P and the force 
moment D = (D x , D z ) is given by |L2], [l3| : 

K _ « 2 (P1-Pft) ^(l-^)(rW + ^) 

L 2 X 3L 2 X ^ 3EL 2 X ' 1 } 

where P\\ ( Pj_) is the component of P parallel (perpendic- 
ular) to the unit- vector n in Fig. |], r is the surface stress 
which generates the z-component of the elastic dipole 
moment (D z = ma), E is the Young's modulus, and v 
is the Poisson's ratio |l7|| . While separating the elastic 
and electrostatic contributions to step interactions is a 
difficult task, we can use the information from Fig. [| to 
obtain estimates for these two effects. 

First, observe that the dipole density due to step- 
rumpling is about 20-30% of the dipole density arising 
from charge spill-out; it is therefore reasonable to cal- 
culate the electrostatic interactions by focusing on the 
latter effect. Next, we see from Fig. || that compared 
to the Sh-step, the Dh-step has a larger value of P\\ due 
to the presence of an additional C atom with dipole p' 
that is aligned closely with the vector n. With increasing 
step- height, the number of such atoms increases, leading 
to even larger values of P\\ and a monotonic decrease in 
the electrostatic contribution (refer to Eq. (||)). However, 
we see from Table. |n| that the step interaction parame- 
ters show a non-monotonic trend, which we attribute to 
an increase in the strength of elastic interactions with the 
height of steps in the following paragraph. 

Using Eq.(|), we calculate the electrostatic contribu- 
tion to the interactions of Sh steps to be ~ leVA, which is 
close to the value of Xf given in Table [□], indicating that 
the step interactions are dominated by the electrostatic 
effects. This observation is also in agreement with the 
fact that a typical value of surface stress r « 100meV/A 2 
yields a small repulsive elastic contribution which is « 
O.OlmeVA. In the case of Dh-steps, a decrease in step in- 
teraction parameter (Ar>) is consistent with a decrease in 
the electrostatic contribution with increasing step- height. 
However, a subsequent increase in the interaction pa- 
rameter (A3) for the Th-step can only be explained if 
the magnitude of the elastic dipole, \D\ = y D 2 + D 2 
increases with the height of the step. In this case, we 
expect both the electrostatic and elastic contributions to 
be comparable in magnitude. 

In summary, by calculating step formation and inter- 
action energies we have shown that Dh and Th steps lead 



to more stable vicinal surfaces, a result consistent with 
experimental observations Jl|-[§. We have described the 
atomic displacements of the surfaces and proposed a sim- 
ple model for the electrostatic interaction of the steps in 
terms of electric dipoles of carbon atoms and step rum- 
pling. 
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